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Abstract 

We use a Glauber-like approach to describe very energetic nucleus-nucleus collisions as 
a sequence of binary nucleon-nucleon collisions. No free parameters are needed: all the 
information comes from simple parametrizations of nucleon-nucleon collision data. Pro- 
duced mesons are assumed not to interact with each other or with the original baryons. 
Comparisons are made to published experimental measurements of baryon rapidity and 
transverse momentum distributions, negative hadron rapidity and transverse momentum 
distributions, average multiplicities of pions, kaons, hyperons, and antihyperons, and zero 
degree energy distributions for sulfur-sulfur collisions at 200 GeV/c per nucleon and for 
lead-lead collisions at 158 GeV/c per nucleon. Good agreement is found except that 
the number of strange particles produced, especially antihyperons, is too small compared 
with experiment. We call this model LEXUS: it is a baseline linear extrapolation of 
ultrarelativistic nucleon-nucleon scattering to heavy ion collisions. 
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1 Introduction 



The program to study the properties of quark-gluon plasma at high energy density is 
in high gear [1J with the construction of the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory (BNL) well underway and due to be completed in 1999. 
In this machine counter-rotating beams of gold nuclei with an energy of 100 GeV per 
nucleon will be collided. The Large Hadron Collider (LHC) at CERN will be completed 
around the year 2005; it will allow for the collision of counter-rotating beams of lead 
nuclei at about 1.5 TeV per nucleon, but will not be dedicated to heavy ion physics as 
RHIC will be. Since 1986 experiments have been performed at CERN's SPS accelerator 
with beams of oxygen and sulfur at 200 GeV/c per nucleon and, lately, lead at 158 GeV/c 
per nucleon, striking fixed targets. In the same time interval similar beams have been 
available at BNL's AGS accelerator at the lower energies of 10 to 14.6 GeV per nucleon. 
It is almost universally accepted that the proper treatment of collisions at RHIC and LHC 
must involve the quark and gluon degrees of freedom. At the AGS hadronic degrees of 
freedom probably suffice (but see 0). The jury is still out concerning collisions at the 
SPS. 

It is oftentimes heard at conferences and workshops that there is a need for a baseline 
calculation of what one would expect at the above heavy ion accelerators if there was no 
new physics; that is, a linear extrapolation of nucleon- nucleon collisions to nucleus- nucleus 
collisions. Construction of such a working model is the goal of this paper. Actually it is not 
so obvious how to make such a linear extrapolation. Nucleon collisions produce mesons, 
and these mesons can collide with other nucleons and mesons, producing an interesting 
cascade of hadrons. We do not consider such a cascade as being a linear extrapolation. 
The collective excitations of such a system are not necessarily trivial, nor is it a simple 
matter to compute or measure all the hadronic cross sections needed to keep track of this 
cascade. Many such models already exist: ARC [1, RQMD @, VENUS §, FRITIOF 
|J, PYTHIA 0, and QGSM being among the most frequently applied to experimental 
data. 

Our interpretation of a linear extrapolation is based on the 40 year old philosophy 
and work of Glauber |J and on the 20 year old rows on rows model of Hiifner and Knoll 
fT0|j as applied to the now disassembled Bevalac at LBNL (Lawrence Berkeley National 
Laboratory). Nucleons from each nucleus follow straight line trajectories, making binary 
collisions with nucleons from the other nucleus. These collisions are as in free space. 
Inelastic collisions produce mesons; the mesons are not allowed to collide with each other 
or with any nucleons. The number of binary collisions suffered by any given nucleon 
depends on the nucleon cross section and on the geometry of the nuclei. The details will 
be given in later sections. It is important to know that this linear extrapolation model 
which we refer to as LEXUS, for Linear Extrapolation of Ultrarelativistic Scattering, has 
no free parameters. 

We will apply LEXUS to published data on S+S and Pb+Pb collisions at the SPS. (At 
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this time the quantity of Pb+Pb data available to us is not as complete as the S+S data.) 
We do not attempt to apply LEXUS to AGS energies. Those energies are probably too 
low to accept the assumption of straightline trajectories as being anywhere near realistic. 

It is important to keep in mind that all LEXUS predictions in this paper are absolutely 
normalized. We have not attempted to tune the results in any sense. 

A conclusion of our paper is that a linear extrapolation of nucleon collisions is con- 
sistent with S+S and Pb+Pb collisions at the SPS in the sense that it gives a good 
representation of the baryon rapidity and transverse momentum distributions, the neg- 
ative hadron rapidity and transverse momentum distributions, the average number of 
pions, and zero degree energy distributions. However, it predicts only about 80% of the 
observed number of charged kaons, 50% of the number of observed neutral kaons and 
lambdas, and 10% of the number of observed antilambdas, all in reference to central S+S 
collisions. This may suggest where new physics lies. 

In section 2 we formulate the model and solve for the basic building blocks, the 
two-particle baryon rapidity distributions. In section 3 we compute the final, observable 
baryon rapidity distribution. In section 4 we compute the baryon transverse momentum 
distribution. In section 5 we compute the average multiplicities of various produced 
hadrons. In section 6 we compute the negative hadron rapidity distribution. In section 
7 we compute the transverse momentum distribution of negative hadrons. In section 8 
we compute the zero degree (calorimeter) energy distribution. Conclusions are drawn in 
section 9. 

2 Formulation and Solution of the Model 

To formulate the model it is convenient to consider a collision between two rows of nu- 
cleons. A nucleus-nucleus collision will be constructed from an ensemble of row-row 
collisions. Refer to the nucleons comprising these rows as projectile and target nucleons. 
Let W mn (yp,yT) represent the 2-particle rapidity distribution for the m'th projectile nu- 
cleon and the n'th target nucleon immediately after their collision. The single-particle 
projectile distribution W^ n (yp) is obtained by integrating the 2-particle distribution over 
the unobserved target rapidity. The index mn then refers to the m'th projectile nucleon 
after colliding with n target nucleons. 



Similarly, W^ n (yr) is the single-particle target distribution obtained by integrating over 
the unobserved projectile rapidity. The index mn then refers to the n'th target nucleon af- 
ter colliding with m projectile nucleons. These single-particle distributions are normalized 
to unity. 




(1) 




(2) 
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Due to the indistinguishability of nucleons the outgoing nucleon with the larger rapidity 
is called a projectile and the other is called a target. 

The 2-particle distribution W mn is obtained by the collision between the m'th pro- 
jectile nucleon, which has suffered n — 1 previous collisions, with the n'th target nucleon, 
which has suffered m — 1 previous collisions. 

W mn (y P , y T ) = J dy'rdy'rW^y^W^MKiy'p + y ' T ^y P + y T ) (3) 

Here we assume that the process is Markovian with kernel K. This is not a necessary 
assumption and could be relaxed. Doing so would result in a correlated cascade. However, 
it would require experimental information on the correlation between the two outgoing 
baryons which is generally not available. 

A basic input is the kernel K which must be taken from experiments on nucleon- 
nucleon collisions. To that end consider a nucleon-nucleon collision in the LAB frame of 
reference so that the initial single-particle distributions are 

W[ (y' P ) = S(y' P -y ) 

WM T ) = %t) (4) 
where j/o is the beam rapidity. Substitution into the evolution equation gives 

W u (y P , y T ) = K(y + -> y P + y T ) . (5) 

Experiments do not measure the correlated two-nucleon distribution over all phase space, 
they only measure the single particle distribution. 

d ^(y) = Wf 1 (y) + W? 1 (y) (6) 

Here the projectile contribution is 

W[ x {y P ) = J dy T W u {yp,y T ), (7) 

and similarly for the target contribution. 

It has long been known that, to good approximation, the distribution of outgoing 
nucleons in a high energy nucleon-nucleon collision is flat in longitudinal momentum or a 
hyperbolic cosine (symmetric about the CM) in rapidity jO, |T^]. This knowledge does not 



uniquely determine the 2-particle kernel K but with the additional, sensible, requirements 
that the projectile distribution be forward peaked and that the simplest mathematical 
representation be used consistent with the data, leads to the parametrization 

K{y' P + y' T ~^yp + yr) = Q(yp - y' T , y'p - yp, y'p - y' T )Q(y'p - vt, yr - y' T , y'p - v't) ( 8 ) 
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where 

Q(s,t,u) = X C -^ + (l-X)5(u). (9) 
sinlit 

In particular the distribution of outgoing projectile nucleons is 

w[M = Q(y, yo, vo-v) = ^^t 1 + A )^° ~ v) • ( 10 ) 

sinhy 



This is the same distribution |Tj| as used in the evolution model proposed by Hwa |L4] to 



describe proton stopping in high energy proton-nucleus scattering and solved and applied 



to data by Csernai and one of the authors The distribution is normalized to unity. 



The parameter A is the fraction of nucleon-nucleon scatterings that result in a hard colli- 
sion and 1 — A is the fraction that are diffractive or elastic. A recent compilation of data 
on p + p — > p + X in the momentum range 12 to 400 GeV/c leads to A = 0.6 []l(J. This 
data and the fit are shown in Figure [IJ. Unless otherwise stated this is the numerical value 
used in the rest of the paper. Physical observables turn out to be rather insensitive to 
small (±0.1) variations in A. There is a small rollover in the data near the projectile and 
target rapidities which is not represented by the parametrization. This has to do with 
precisely how one separates hard inelastic and diffractive collisions. Our results are only 
as good as the input parametrizations; in the future it might be worthwhile to treat these 
components on a finer level. 

There is an obvious and useful symmetry between the single-particle projectile and 
target distributions. 

wL(y) = wLivo - y) (n) 

These distributions are not independent. In the present formulation of the model only 
single particle observables may be reliably computed. Hence we only need to compute 
the W^ n (y). For this we need an evolution equation. It is obtained by integrating eq. (|3|) 
over the target rapidity and using eq. (Bf) for the kernel. 



W. 



inn 



(y) = J dypdyTW^ n _ l {y P )w^ l {y - y T )Q{y - yr, yp -yr,y- vp) (12) 



It only remains to solve this Boltzmann-like equation. This can be accomplished by 
starting with the initial distribution eq. (|]) and then iterating over all m and n. 

A closed form expression can be given for the first nucleon in the row undergoing an 
arbitrary number of collisions [15]. 

-I k-l 



smh y 




In 



/ sinh y 
\ sinh y 



+ (l-A)%-») (13) 



Unfortunately we were not able to find a closed expression for arbitrary mn, so we solved 
the equations numerically for m and n up to and including 14. The solutions have the 
form 

WLiv) = W P mn {y) + (1 - X) n 6(y - y) (14) 
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where W is a continuous function albeit with logarithmic singularities at y = and 
y = yo- These very soft singularities are a consequence of the explicit function Q chosen 
above. When comparing with experiment it should be remembered to smooth these by the 
experimental resolution of the detectors. See also Figure [I] and the previous discussion. 
We have constructed numerical tables of W. These may be obtained from the web site of 



one of the authors [17]. Some representative examples are plotted in Figure 



3 Baryon Rapidity Distribution 

In this section we will apply the most obvious output of the model, the baryon rapidity 
distribution, to the available experimental data on nucleus-nucleus collisions. But first we 
must describe how to make a nucleus-nucleus collision out of row on row collisions. This 
is standard material for any Glauber-like model so we shall go over it without too much 
discussion. 

Consider a collision between a projectile nucleus and a target nucleus with an impact 
parameter b. We can think of this approximately as a sum of independent collisions 
between rows of nucleons as illustrated in Figure |3|. Two rows will collide when the 
transverse position Sp of the projectile row relative to the projectile nucleus' center of 
mass and the transverse position of the target row relative to the target nucleus' 
center of mass are related by sp = b + sp. The average number of nucleons in each row, 
with cross sectional area equal to the nucleon-nucleon cross section, are 

z/p(sp) = a NN J dzp P (s P ,z) 

u t (st) = a NN J dzp T {s T} z) (15) 

where p is the baryon density and z is the longitudinal coordinate. Let V?(s P ) and 
V?(s T ) denote the probability of having fh and n nucleons in the projectile and target 
rows, respectively, at the given impact parameter. We defer the actual choice of these 
probabilities. 

There will be fluctuations in the number of nucleons in each row. Similarly, there will 
be fluctuations in the number of collisions suffered by any given nucleon in a row. Taking 
these into account results in the contribution of the projectile nucleons to the final baryon 
rapidity distribution as follows. 

rIM A P rh A T i2 „ 

>b)=EEE wUv) / —vlMvgM (ie) 

ay M=lm=ln=Q J (Jnn 

There is an analogous expression for the target contribution. The total rapidity distribu- 
tion at the fixed impact parameter is the sum of the projectile and target contributions. 

Now the choice of the Vs must be made. One candidate is a Poisson distribution. 
Its disadvantage is two-fold: it overestimates the magnitude of the fluctuations (they are 
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unlimited in a Poisson but are limited in reality by the number of available nucleons) and 
it would require knowledge of the W£ n for arbitrarily large values of m and n, which is 
computationally infeasible. A natural alternative is a binomial; but with what maximum 
value? It is physically unreasonable to allow all nucleons in a nucleus to fluctuate into one 
row. The maximum average value of the nucleons in a given row is about 10 in a large 
nucleus. We have chosen the maximum value to be N max to be 14. Varying this number 
by 1 or 2 does not change any distribution by more than a percent. Explicitly we have 
chosen 



KM 



( A^max \ 


v P (sp) 


m 


1 - 


U P (Sp) 


\ fn J 


N 

1 v max 






A^max 



6(iV max -m). (17) 



Here is the step function. For the nucleon-nucleon cross section we use a constant value 
of 4 fm 2 which is appropriate for beam energies ranging from tens to hundreds of GeV. 
For the nuclear density distribution in lead we use a Woods-Saxon 



p(r) 



1 + exp[(r — b)/a] 



with parameters a = 0.546 fm and b = 6.62 fm. Normalization to A = 208 fixes po 
0.1604 fm -3 . For sulfur we use a three-parameter Gaussian 



p(r) = ^ o(1 + Wr2/&2) (19) 
PK ' l + exp[(r 2 -6 2 )/a 2 ] K ' 

with parameters w = 0.160, a = 2.191 fm, and b = 2.54 fm. Normalization to A = 32 
fixes po = 0.226 fm -3 . 

To compare with experimental measurements we must know the trigger conditions; 
that is, we must know with what probability any particular impact parameter is accepted 
by the detector. Such a trigger can best be accommodated by an event generator, for then 
the output of the theory can be sent through the experimental filter. Since our model is 
not an event generator, at least in its present form, we can only attempt to simulate the 
trigger as best we can. The procedure we shall use is to make a sharp impact parameter 
cutoff. The total nucleus-nucleus cross section may be computed following Karol |[L9||. 



tot 
a A P A T 



J d 2 6{l-exp[-/(b)]} (20) 
Here / is the geometrical overlap function of the two nuclei. 

/( b ) = / ^^ p (spK(s T ) (21) 

J &NN 

If only those nucleus-nucleus collisions are accepted with an impact parameter less than 
or equal to b cut then the corresponding cross section is 



<?A P A T (b C ut) = J d 2 b{\ - exp [-/(b)]} 6(6^ - b) . 



'.22) 
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Clearly 

a% AT = lim a ApAT (b cut ). (23) 

The impact parameter cutoff can be adjusted to reproduce a given centrality cut. For 
example, if only the 6% "most central" collisions are accepted then b cut is adjusted such 
that o ApAt (6 cut ) = 0.06(7^* At . 

Now we compare with experiment. Figure ^ shows the proton rapidity distribution 
measured by NA35 [^(J for the 2% most central collisions of S+S at a beam energy of 
200 GeV/c per nucleon. The solid symbols are the actual measurements and the open 
symbols are those obtained by reflection about midrapidity. Except possibly for the 
nuclear fragmentation regions, that is, within ±1 unit of rapidity of beam and target, the 
data is represented very well by LEXUS. 

The first measurements of the proton rapidity distribution in Pb+Pb collisions at 158 



GeV/c per nucleon have just recently been published by NA44 They are shown in 



Figure | and represent the 6.4% most central collisions. The systematic plus statistical 
uncertainties together are very large, and there are only two measured points. To compare 
with the proton distributions in Pb+Pb we should point out an uncertainty in our current 
application of LEXUS. We have not distinguished between outgoing protons and neutrons 
but have only counted baryons. For collisions of charge asymmetric nuclei at high energy 
it is reasonable to expect that nearly 1/2 of the outgoing baryons will be protons due to 
the preference to convert more neutrons into protons than vice versa. This simply follows 
from phase space and entropy. Accepting that, the LEXUS predictions are in very good 
agreement with the data. 

So far NA49 has not published its measurements of the baryon rapidity distribution in 
Pb+Pb collisions. However, we would like to point out a curious feature of the symmetry 
between forward and backward going particles in the center of mass frame. In an ideal 
measurement it would be symmetric in collisions of identical nuclei when averaged over 
a sufficient number of collisions. But if there is a zero degree calorimeter used to select 
on central collisions then any particles going into that calorimeter may be considered 
as lost and not counted in the rapidity distribution. This induces a forward/backward 
asymmetry. The effect is shown in LEXUS in Figure |] (details will be discussed in section 
8). The upper symmetric curve represents the proton distribution without a zero degree 
calorimeter. The solid curve presents the depletion of protons in the forward direction. 
It is clear that different results are obtained if one measures only for y > 3 and reflects 
about y — 3 than if one measures only for y < 3 and then reflects. 



4 Baryon Transverse Momentum Distribution 

Even though LEXUS assumes straightline trajectories it is still possible to get an enhance- 
ment of the baryon transverse momentum compared to nucleon-nucleon collisions. Every 
time a collision occurs the baryons can get a transverse kick. As long as they continue to 
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travel with high velocity this will not negate the assumption of straightline trajectories. 
The path of the baryon may be thought of as a random walk in transverse momentum 
space. The average transverse momentum-squared of a baryon is related to the average 
number k of collisions it has suffered according to 

(Pr)* = k(p 2 T ) NN (24) 

where (p\) nn is the average in a nucleon-nucleon collision. This quantity is approximately 



beam energy independent for the beam energies of interest to us |]22| , |23|| . We will use 
(p 2 T ) NN = 0.282 (GeV/c) 2 . 

We will assume that the rapidity and transverse momentum distributions factorize in 
the elementary nucleon-nucleon collisions, which is close to fact except near the edges of 
phase space. The transverse momentum distribution is well represented by the thermal 
form 

= MnwrK^mv/Tx) (25) 

where my = J m 2 N + p\ is the proton transverse mass and M\ normalizes the distribution 
to one. Taking {p\)nn = 0.282 (GeV/c) 2 converts into T\ = 113 MeV. After suffering k 
collisions the baryon transverse momentum distribution becomes broader as reflected in 



a higher temperature Tj, [24] which is determined by 



rod 

M k j dp T p%m T Kx(m T /T k ) = k(p^) NN . (26) 
Jo 

The transverse momentum distribution of projectile baryons can now be computed by 
averaging over the number of collisions. For y ^ y , p? ^ it is 

x / d \ M b M -b) f^VZ( ST m SP ). (27) 

J ° A P A T {°cut) J &NN 

The reason for the sum over k at fixed n is that the probability of any given scattering 
to be considered hard is A, which is binomially distributed. The k = term just results 
in pt = 0. There is an analogous expression for target nucleons. 

Now we compare with experiment. The transverse momentum distribution in the 2% 
most central S+S collisions at 200 GeV/c per nucleon has been measured by NA35 



The data span the rapidity range 0.2 < y < 3.0 and are shown in Figure 0. The shape 
and the absolute normalization are in very good agreement with the results of LEXUS. 

The transverse mass distribution for protons from the 6.4% most central Pb+Pb 
collisions at 158 GeV/ c per nucleon have been measured by NA44 JH]]. The measurements 
were performed at y = 2.10 and at y — 2.65. They are shown in comparison with 
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LEXUS in Figure ^. The higher rapidity data is again in very good agreement with 
LEXUS; the lower rapidity data is in less good agreement. The most likely explanation 
is that we have assumed that elastically or diffractively scattered nucleons acquire zero 
transverse momentum. Allowing for it would increase the distribution at small px nearer 
the projectile and target rapidities. This is a topic for future investigation. 

5 Average Multiplicities of Produced Hadrons 

The production of secondary hadrons, such as pions and kaons, or the conversion of the 
incident nucleons to hyperons, carries additional information about the collision dynamics, 
in particular the entropy. In LEXUS we assume that mesons are created in the collisions 
of the cascading baryons as in free space. Once created, the mesons are assumed to freely 
disperse: they no longer participate in the nuclear collision. The incident nucleons may 
very well change their character as they cascade: protons may convert to neutrons and 
vice versa, or they may be excited into various A, iV*, or hyperon states. In this section 
we will compute the average multiplicities of: charged hadrons h~; kaons K + , K~ , Kg] 
hyperons A and A. We could compute the full multiplicity distributions, but in this paper 
we are content to get the average multiplicities. 

Let {X(s))nn represent the average number of mesons of type X produced in a 
nucleon-nucleon collision at center-of-mass energy ^fs. The average number of such 
mesons produced in a collision between the m'th projectile nucleon and the n'th tar- 
get nucleon is 

(X mn ) = \J dypdy T W^ n _ 1 (yp)W^_ ln (y T )(X(yp - y T ))NN , (28) 

where = 2m n cosh [(yp — yp) / 2]. The total number produced in a given row-row 
collision is obtained by summing over all nucleon-nucleon collisions. 

Ap At fh n 

(x(s P ,s T )) = EE^l^M E E(^-) (29) 

m=l n=l m=l n=l 

Finally we need to sum over all rows and over all allowed impact parameters. 

(X(b < b cut )) = [ 7j-^®(b cut - b) [ -^{X(bp, s t )> (30) 

J ° ApA T {Ocut) J &NN 

This can be written as a trace over the product of two matrices. 

Let us first consider the production of negatively charged hadrons. Gazdzicki has 
shown []3SJ that the average multiplicity of charged hadrons in isospin averaged nucleon- 
nucleon collisions for laboratory momenta ranging from 2 to 400 GeV/c is fit to within 
about 6% by the simple parametrization 

(h-) NN = 0.78AF NN7T (s) (31) 
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where -Fawv is the Fermi variable modified for the pion production threshold 



(y/s - 2m N - m n ) 3/4 

'AWttISJ = 773 i 62 ) 



and where y/s is given in GeV. There are two caveats to using this parametrization in 
LEXUS. The first is that for light nuclei such as oxygen or sulfur the isospin averaging 
is alright, but for heavier nuclei such as gold or lead it is not. However, after one or two 
collisions protons are more likely to transform into neutrons than vice versa because of 
phase space or entropy. Then isospin averaging becomes a better approximation. The 
second is that occasionally a nucleon can convert into a hyperon. There is essentially 
no experimental information on the charged hadron multiplicity in a hyperon-nucleon 
collision, for good reason. Therefore, for lack of any other information, we shall continue 
to use Gazdzicki's parametrization. 

There is somewhat less experimental information on the multiplicity of kaons pro- 
duced in nucleon-nucleon collisions. These data have been compiled by Gazdzicki and 
Rohrich |26| and shown to be rather well-defined functions of the Fermi variable. The 



parametrizations we shall adopt are given below. 

(K + ) NN = 
(K-) NN - 



n P 2 

NNK+ 



b+ (F NNK + - c) 2 
nF 3 

ar NNK- 

b + {F NNK - - c) 2 

^ = » + (w,> (33) 

The parameters a, b, c are displayed in Table 1. 

The experimental data on A and A production was also analyzed by Gazdzicki and 
Rohrich. We have constructed the following parametrizations for use in LEXUS. 



(A) 



NN 



n 

NNA 



b + (F NNA - cf 



These parameters a, b, c are also listed in Table 1. 
In all cases above the F^nx are defined as 

p iVl - M x y/± 

" NNX — -^Jl , [OO) 

where the threshold energies are: Mx = + M\ + for X = A, K + and K§; 

M x = 2M N + M K + + M K - for K~\ and M x = 2M N + 2M A for A. We have taken 
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into account the uncertainties in the production cross sections in the elementary nucleon- 
nucleon collisions as illustrated in Figure |9| for the A hyperon. In addition to making a 
best fit to the data, upper and lower envelopes are constructed which roughly pass through 
the upper and lower error bars, respectively. (The parameters for these envelopes are not 
tabulated here.) 

Calculated results for S+S collisions are shown in Table 2 and compared with mea- 
surements of NA35 [27L |28] for 2% centrality. The LEXUS results are given for 0%, 2% 



and 4% most central collisions based on impact parameter. These represent the variation 
in computed abundances with small variations in impact parameter away from 6 = 0. The 
experimental data is based on a definition of centrality as those collisions with the highest 
value of transverse energy. This does not exactly correspond to a sharp range of impact 
parameter because of fluctuations. A better comparison would involve impact parameter 
smearing. That is outside the scope of this paper. 

The average number of negatively charged hadrons predicted by LEXUS is in very 
good agreement with the measurements. The average number of positive or negative kaons 
obtained from LEXUS is about 80% of that observed. This is a two standard deviation 
effect when account is made of the uncertainty in the production rate in nucleon-nucleon 
collisions. The number of short-lived neutral kaons is a factor of 2 too small, as is the 
number of lambdas. The number of antilambdas is nearly an order of magnitude too 
small compared to experiment. This almost certainly indicates a failure of the model, 
most likely as a result of the neglect of multiple scattering of produced mesons. 

Table 3 is a prediction of LEXUS for the 5% most central collisions of Pb+Pb at 158 
GeV/c per nucleon. No data has yet been published for these abundances. 



6 Rapidity Distribution of Secondaries 

It is important to know where the produced particles emerge in rapidity space. Several 
types of detectors are able to measure the negative charged hadron rapidity distribution 
d(h-)/dy. 

In nucleon-nucleon collisions the charged particle rapidity distribution is approxi- 
mately Gaussian 



d(h )nn (h )nn 

exp 



dy V27TO-L(Z/rel 

with a width given by the Landau model 



-(y-ycm) 2 / 2cr i(2/rel) (36) 



, iforel)= 8_j_, ln (_g (37 ) 

where y re \ = yp — yr is the relative rapidity of projectile and target and y cm = | (yp + yx) 
is the rapidity of the center-of-mass. Cq is the speed of sound of the produced matter. 
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Not surprisingly Cq ~ 1/3 provides a good fit to nucleon- nucleoli data P9"| . We have used 
a very slightly smaller value of 0.32 corresponding to a free gas of massive pions at a 
temperature of 140 to 160 MeV. The distinction between 1/3 and 0.32 is irrelevant for 
the purposes of this paper. 

It is straightforward to compute the distribution d(h~)/dy in nucleus-nucleus col- 
lisions. All we need to do is replace the quantity (X) NN with expression flHEP in eqs. 



The rapidity distribution of h has been measured in S+S collisions by NA35 



The results for the 2% most central collisions are shown in Figure |10| The prediction of 
LEXUS is also shown and presents a very good representation of the data. Predictions 
for the 5% most central collisions of Pb+Pb at 158 GeV/c per nucleon are shown in 



Figure 11. No data have yet been published. 



7 Negative Hadron Transverse Momentum Distribu- 
tion 

Mesons are produced when nucleons collide, and since the colliding nucleons are under- 
going a random walk in transverse momentum, the mesons will acquire extra transverse 
momentum too. This is taken into account in the following way. 

Suppose that the center-of-mass frame takes discrete steps v in transverse velocity 
when baryons undergo hard scattering. The magnitude is fixed but its direction in the 
transverse plane is random. Then one can show that the average transverse momen- 
tum squared of pions produced in a collision between two nucleons which have together 
undergone i previous scatterings is 



where 



Mi = (Pt)tt + 



1 -v 2 



\(p 2 T ) n + ml 



i-l 

E 

k=0 



(3* 



(39) 



Here (p 2 -)^ is the average squared transverse momentum of pions in elementary nucleon- 
nucleon collisions. Identifying 

AmW (40) 



(Pt)nn 

i — v~ 

we obtain v = 0.272. 

The average value of the pion transverse momentum in elementary nucleon-nucleon 
collisions varies surprisingly little for beam momenta ranging from 11.6 GeV/c to 195 
GeV/c |22|| . We shall use (pt)^ = 0.155 (GeV/c) 2 . For pions, too, a good representation 
of the transverse momentum distribution in nucleon-nucleon collisions is the thermal form 
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(P5|) with the nucleon mass replaced by the pion mass. This results in a pion temperature 
in nucleon-nucleon collisions of 133 MeV. 

After % previous scatterings of nucleons the pion distribution becomes broader with a 
temperature Tj determined by 

Mi dp TP 3 T m T K 1 {m T /T i ) = {p 2 T ) i . (41) 
Jo 

The full distribution of pions in rapidity and transverse momentum is then computed 
analogously to the procedure in section 5. 



A P A 



p T dp T dy 



J A.J-' J i-± lib I ti n 

E E f^MKM E E A / dypdyrW^MW^M 

— 1 = 1 „ 1 1 J 



X 



m=l n=l m=l n=l 

(n-(y vc i)) NN exp - (y - y cm ) 2 /2a 2 L (y Tel ) 



V2na L (y re i) 

x M n+m -2m T K l (m T /T n+m _ 2 ) . (42) 



A similar analysis can be done for kaons. In this case we use p2| (pt)k — 0.290 (GeV/c 



\2 



We assume that the sum of the K~ and ir~ multiplicities equals the h~ multiplicity. 

In Figure [L2 we show the h~ transverse momentum distribution for the 2% most 
central S+S collisions for the rapidity interval 0.8 < y < 2.0. The data is from NA35 
|20[| . In Figure [13| we show the distribution for the interval 2 < y < 3. The agreement 



is acceptable with no surprises. The slight underestimate at very small pt could be a 
result of too crude an approximation to the transverse momentum distribution in nucleon- 
nucleon collisions. It could also be a result of multiple scattering among the produced 
pions ]3Tfl. We do not attempt to compute the distribution for p T > 1.5 GeV/c since 



the parametrization chosen is not representative of high p? pion data in nucleon-nucleon 
collisions. 



8 Zero Degree Energy Distribution 

Many experiments have what is known as a Zero Degree Calorimeter (ZDC) which mea- 
sures the energy carried by forward going particles in a collision. It is in some sense a 
measure of impact parameter since any projectile nucleons which participate in the colli- 
sion get scattered away from the forward direction; to first approximation only spectator 
nucleons go into the ZDC. There is a monotonic relationship between the average number 
of spectator nucleons and the impact parameter. The ZDC is oftentimes used as a cen- 
trality trigger whereby acceptance of only small energy deposition roughly corresponds to 
central collisions. 

The physics of the ZDC is much more complicated than the basic idea presented 
above. The ZDC has a very specific response to a given hadron (p, n, it , K~, etc.) 
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with a given LAB momentum. This response must be accounted for to get an accurate 
comparison between a model calculation and the data. In this paper we shall make a zero 
order estimate of the ZDC energy distribution for heavy ion collisions, and then make a 
first order correction. 

In a heavy ion collision a certain number of projectile nucleons will not scatter but 
will continue along the beam direction without deflection. Knowing this number we can 
compute the energy deposited in the ZDC (see caveat above). At a given impact parameter 
the average number of projectile spectator nucleons is 



N p , 

- 1 * spec v 



(b) = I — E*)m^(Br)(l- A)». (43) 
J cr NN mn 



Here rh is the number of projectile nucleons in the row with probability distribution 
and V^(l — A) n is the probability of encountering n target nucleons without suffering any 
hard scattering. All possibilities are summed over as are all rows. With the probability 
distribution (|I7|) we get 



•max 

* max I 



< ec (6) = f ^p(sp) [1 - Xp t (s t )/N b 
J (Tatty 

~ [— -vp(sp)exp[-\v T (s T )} . (44) 

J &NN 

This has the interpretation of the average number of projectile nucleons in a row times 
the probability of not making a hard collision with any of the target nucleons, integrated 
over all rows. When A — > or when b — > oo the average number of spectator projectile 
nucleons approaches Ap as it should. 

There will be fluctuations in the number of spectator nucleons even at fixed impact 
parameter. The dispersion can be computed in the same way as the average number. 

d 2 sp 



D 2 spec (b) = [ ^Y.KMmV T n {s T ){l-\Y[l-{l-\Y] 

J 0" NN ™„ 



&NN mn 
12 



^z/p(sp) {[1 - Az/ T (s T )/iV max ] 7V — - [1 - A(2 - A)z/ r (s r )/iV max ] JV — } 
a NN <• 

~ / ^Msp) {exp [-XvtM] - exp [-A(2 - A)z/ r (s r )]} (45) 

When A — > or when b — ► oo the dispersion goes to zero, consistent with every single one 
of the projectile nucleons entering the ZDC. 

We will use the central limit theorem to approximate the conditional probability 
distribution of the energy carried by projectile spectators into the ZDC at fixed impact 
parameter by a Gaussian. 

;t^(£zdc|&) = * ex P {" [£zdc - (E(b))} 2 /2D 2 (b)El am } (46) 

dE ZDC V27rD(b)E he£Lm 1 J 
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Here we have made the identification (E(b)} = E^e^N^^b) and similarly for the disper- 
sion. To get the ZDC energy distribution we integrate over impact parameter. 

f d 2 bJ^(E ZDC \b) (47) 



dEzBC J dEzBC 



Now we compare with experiment. Data from NA35 for S+S |3l] is shown in Figure [14 
as well as the the cross section for forwardgoing nucleons in LEXUS (dashed curve). The 
central plateau is just about right, as it should be, reflecting the basic geometry of the 
nuclei. The data go beyond the kinematical limit of 6.4 TeV in a single S+S collision 
as a result of inefficiencies in the detector. LEXUS predicts too much cross section for 
forwardgoing energies of 1 to 2 TeV. 



Data from NA49 for Pb+Pb collisions [32 is shown in Figure 15. Once again the 



central plateau is correctly reproduced, but the shoulder at low energy, indicating the 
most central collisions, is shifted too far left by about 4 TeV. This ZDC is improved over 
that of NA35 and does not go beyond the kinematic limit of 32.9 TeV. 

In actual experiments there is a finite opening angle 8q for particle acceptance in the 
ZDC. This is generally a fraction of a degree, in the LAB frame of course. Some of the 
hard scattered nucleons, both projectile and target, may emerge with a LAB angle smaller 
than this. This effect will tend to increase the energy flow into the ZDC at a fixed impact 
parameter. The additional energy entering the ZDC is 

r d 2 N p 
E <e<6 (b) = J dydpr (m T cosh y - m N ) ^ (b, p T , y) 

x (# — tan -1 (p T /?T7 T sinhy)^ (48) 

Here d 2 N p /dpxdy is the same as expression (p7f) but without the averaging over impact 
parameter. Then in eq. fllB]) we identify (E(b)) = i?beam-V S pec(^) + -^o<e<e (^)- 111 addition, 
the dispersion in the energy entering the ZDC increases and can be estimated by 



D 2 {b) = fflL 5>„(s p ) V n (s T ) £ a 2 mn . (49) 

J &NN ™~ 1 



mn m=l 



Here 



2 

° mn 



dy / dp T E 2 K ^-^(pr, y) - dy dp T E K -^(pt, y) 
Jo dprdy \J Jo dpxdy J 

+ Target Contribution (50) 
where Ejc = (mj'cosh.y — m^), p™ ax = tan^m^ sinhy, and 

+ (l-X)-8{y-y )8(p T ). (51) 
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The target contribution is obtained by simply changing y — > yo — y. 

Now we recompare with experiments. The results of LEXUS with an opening angle 
of 0.15° is shown in Figure [14] by the solid curve. The result of allowing some of the 
scattered nucleons to enter the ZDC is obvious; it reduces the cross section at small 
forward energies. The results for Pb+Pb with an opening angle of 0.3° are shown in 
Figure also with a solid curve. The better agreement for Pb than for S is somewhat 
mysterious. The opening in the ZDC for both experiments was actually a square, not a 
circle, with an aperture of 86 //sr. It is quite possible that the simple angle cutoff we have 
used does not do justice to the complicated workings of these calorimeters. It is also quite 
possible that our use of a baryon distribution which factorizes in rapidity and transverse 
momentum in nucleon-nucleon collisions is too crude near the edges of phase space, such 
as for very forwardgoing nucleons. 



9 Conclusion 

In this paper we have constructed a means to make a linear extrapolation of nucleon- 
nucleon collisions to very high energy nucleus-nucleus collisions. We call this extrapolation 
procedure LEXUS. There is no reference to quarks, gluons, strings, Pomerons, or QCD. 
The treatment is simply based on a sequence of binary nucleon-nucleon collisions as in 
free space. We know that this treatment cannot be exact as all of our accumulated 
knowledge of QCD and high energy physics of the last twenty-five years can attest. But it 
is important to do these calculations as a baseline against more detailed, but of necessity 
approximate, treatments based on perturbative and nonperturbative QCD to discover 
thermodynamic properties of quark-gluon plasma and hadronic matter. What have we 
learned? 

The rapidity and transverse momentum distributions of baryons in central sulfur- 
sulfur and lead- lead collisions at the SPS seem to be well-described by LEXUS. The same 
can be said for the multiplicity of negatively charged hadrons and their rapidity and 
transverse momentum distributions. The zero degree energy distribution for lead-lead 
comes out just about right; for sulfur-sulfur the agreement is less good. This may be due 
to our rather simple treatment of elastic and diffractive nucleon-nucleon collisions which 
are more likely to influence the outcome of the collisions between smaller nuclei. 

What is not reproduced so well by LEXUS is the abundance of strange particles. 
Formation of quark-gluon plasma increases the number of strange quarks in the system 
due to its relatively small mass in comparison to the kaon mass [j33fl . But it is also 



possible to increase the number of kaons over that produced in nucleon-nucleon collisions 
by multiple scattering and attendant pair production of kaons [|3~4" |. Most notably the 



number of antilambdas is too small in LEXUS in comparison to central sulfur-sulfur 
data by almost an order of magnitude, indicating that there are proportionately many 
more antistrange quarks in the heavy ion collision than in nucleon collisions p5| . This is 
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interesting physics. It has been discussed several times at Quark Matter conferences |TJ. 
Our calculations confirm it. We anxiously await published data on strangeness in central 
lead-lead collisions. 

There is a preliminary conference report on a forward - backward azimuthal asymme- 
try in Pb-Pb collisions |36] . This asymmetry is similar to that seen at much lower Bevalac 
energies. It has been interpreted variously as a collective bounce-off, as if the two nuclei 
were behaving as fluids, and as the absorption of particles in the cold spectator matter. 
In either way of thinking the present linear extrapolation of nucleon-nucleon scattering 
does not take this crosstalk of rows into consideration. It would be challenging to do so 
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A linear extrapolation like LEXUS is only as good as the data input from nucleon- 
nucleon collisions. In this paper we have used reasonably accurate yet simple parametriza- 
tions of the basic input data. Nevertheless improvements can be made. For example, our 
treatment of the elastic and diffractive components of nucleon-nucleon collisions could be 
represented more accurately but at the cost of significant complication to the solution 
to the model as described in section 2. A better prediction for the baryon momentum 
distribution near the projectile and target rapidities would likely result, as would a de- 
scription of the energy deposited in a zero degree calorimeter. Even then, our ability to 
make a linear extrapolation will be hindered by the lack of experimental measurements 
of many-particle correlations, such as between the two outgoing baryons in an elementary 
collision, or the correlation between the rapidities of outgoing baryons with the number of 
produced mesons. It is unlikely that all the exclusive cross sections for nucleon-nucleon, 
nucleon-hyperon, and hyperon-hyperon collisions at all the energies of relevance will ever 
be known experimentally. This is a particular shortcoming if one wants to make a Monte 
Carlo event generator out of LEXUS. In case of improvements to LEXUS in the future 
the model described in this paper will be known as version 1.0. 

We have not made any comparison to transverse energy distributions as measured 
by electromagnetic calorimeters. Many such measurements have been made, but it is a 
challenging task to normalize them in such a way that we can know how much energy 
each type of hadron (p, n, tc°, K~ , etc.) actually deposits to within a certain accuracy 
to make a comparison meaningful. This will be the subject of a future application of 
LEXUS. 

Also under study is an application of LEXUS to some very interesting results on 
the production of J /if), photons, and lepton pairs as discussed in recent Quark Matter 
conferences JTJ. 

Finally, it would be interesting to make LEXUS predictions for RHIC where we cer- 



tainly anticipate it to fail, lacking a description of hard perturbative QCD [38|, |39| and the 



inside-outside cascade effect ||40|| . The difficulty is that there is no nucleon-nucleon data 
available at RHIC energy before the turn-on of the accelerator. Nevertheless a prediction 
for the baryon rapidity distribution may be considered and is now in progress. 
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0.00203 


b 


6.92 


2.62 


4.28 


4.03 


5.67 
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2.24 


2.20 


0.794 


1.16 


2.14 



Table 1: Table of the parameters for the fits in eqs. (|33|-p4]) 
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Table 2: The average particle multiplicities in S±S collisions. The NA35 data should be 
compared against the 2% most central collisions in LEXUS. 
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Table 3: Predictions for 158 GeV/c Pb±Pb collisions assuming and 5% centrality. 
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Figure 1: Proton rapidity distribution in p + p — > p + X in the center of mass frame. From 
the top: pi a b = 12,24,69, and 400 GeV/c. The filled symbols indicate experimentally 
measured data points and the empty symbols are the reflected ones. The solid lines are 
dN/dy = Acosh(y)/ sinh(?/o/2) with corresponding maximum rapidity y . The data were 
assembled in Ref. [[UJ. 
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Figure 3: A schematic view of a row on row collision. 
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Figure 4: Proton rapidity distribution for a S+S collision at 200 GeV/c with a centrality 
of 2%. The solid line represents our calculation. Filled diamonds are data from NA35. 
Open diamonds and the dashed line are the reflection of the left half. 
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Figure 5: Proton rapidity distribution for a Pb+Pb collision at 158 GeV/c with a central- 
ity of 6.4%. The solid line is our result. The data points from NA44 are marked by filled 
diamonds. The open diamonds and the dashed lines are the reflection of the left half. 
Error bars on the NA44 data represent statistical errors, the limiting short bars represent 
systematic errors. 
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Figure 6: Proton rapidity distribution for a Pb+Pb collision at 158 GeV/c with 5% 
centrality. The solid line is our result, the dashed line is the reflection about the mid- 
rapidity. The asymmetry is caused by the high rapidity particle loss to the zero degree 
calorimeter. 
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Figure 7: Proton transverse momentum distribution for a S+S collision at 200 GeV/c 
with a centrality of 2%. The solid line represents our calculation. Data are from NA35. 
Rapidity range is 0.2 < y < 3.0. 
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Figure 8: Proton transverse momentum distribution for a Pb+Pb collision at 158 GeV/c 
with a centrality of 6.4%. Filled diamonds are NA44 data and the solid curve is our result. 
The upper curve is for y = 2.65. The lower curve is for y = 2.10 scaled by a factor of 
1/10. 
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Figure 9: Data and fits for A hyperon production in p + p collisions. The middle line is 
the weighted fit for the whole data, the upper and the lower lines represent maximum and 
the minimum envelopes. Experimental data are taken from Ref. . 



31 




y 

Figure 10: Rapidity distribution of h~ in 200 GeV/c S+S collisions with a centrality of 
2%. The solid line represents our result. Filled diamonds are data from NA35. Open 
diamonds and the dashed line are the reflection of the left half. 
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Figure 11: LEXUS prediction for the rapidity distribution of negatively charged hadrons 
in 158 GeV/c Pb+Pb collisions with a centrality of 5%. 
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Figure 12: Negative hadron transverse momentum distribution for a S+S collision at 200 
GeV/c with a centrality of 2%. The solid line represents our result and the filled diamonds 
are NA35 data. Rapidity range is 0.8 < y < 2.0. 
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Figure 14: The zero degree energy distribution for 200 GeV/c S+S collisions. The solid 
line represents LEXUS with an opening angle of 0.15 degree. The dashed line represents 



LEXUS with only the spectator nucleons. Data are from NA35 J3TJ . 
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Figure 15: The zero degree energy distribution for 158 GeV/c Pb+Pb collisions. The solid 
line represents LEXUS with an opening angle of 0.3 degree. The dashed line represents 



LEXUS with only the spectator nucleons. Data are from NA49 |32 
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